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Abstract 

We extract the running coupling based on Creutz ratios in SU(2) lattice gauge theory with two 
Dirac fermions in the adjoint representation. Depending on how the extrapolation to zero fermion 
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coupling. This behavior is consistent with other findings which indicate that this theory has an 
infrared fixed point. 
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I. MOTIVATION 



Understanding the allowed model space in extensions to Standard Model is a research 
program that gains motivation at the onset of the Large Hadron Collider. In technicolor 
models, the Higgs mechanism occurs through condensation of new fermions that are subject 



to a gauge interaction that is strong at the TeV scale jl,|2|. Walking technicolor is a version 
of this theory that can address certain difficulties associated with flavor-changing neutral 
currents js-Oj. Higher representations of the gauge group are believed to avoid problems with 



the S-parameter, i.e. electroweak precision constraints |10l. All of this has motivated the 
study of Minimal Walking Technicolor 12| using lattice techniques, in order to see where 
this theory lies with respect to the conformal window. In other words, there is no guarantee 
that the theory actually walks; finding out one way or another is to solve a strong coupling 
problem. Here we discuss a method, based on the behavior of Creutz ratios, that does just 
that. The approach taken here is different from the Schrodinger functional method that 
was employed for SU(3) gauge group with triplet fermions i n |l3| . SU(3) gauge group with 



sextet fermions [14!], and Minimal Walking Technicolor [15-17|. The latter works indicate 
the existence of an infrared fixed point. We obtain results consistent with that conclusion 
by this alternate Creutz ratio approach. 



II. THE METHOD 



In 



18| a method for measuring the running gauge coupling on the lattice using Creutz 



ratios was proposed. In [l9|, |20| this method was apphed to pure SU(3) lattice Yang-Mills, 



and in 



20 



2l| this method was applied in a pre liminary way to SU(3) gauge group with 



sixteen triplet fermions, and more recently in [22j to SU(2) gauge group with eight fiavors 
of fundamental representation fermions. We will use this approach to search for the phe- 
nomenon of backward running, which is a smoking gun for a nontrivial infrared fixed point, 
since forward running is assured at weak coupling by perturbative methods. We apply the 
method of Creutz ratios to the theory of Minimal Walking Technicolor: SU(2) gauge the- 
ory with two fiavors of fermions in the adjoint representation. Our lattice action is Wilson 
fermions and a plaquette gauge action. 
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A. Outline 

The proposal [18|] defines the running coupling g{L) at a scale L = Na associated with 
the spatial extent of the lattice, where N is the number of sites in spatial directions and a 



is the lattice spacing.^ This is done through the Creutz ratio [23|: 



where W{I, J) is the expectation value of the trace of the rectangular I x J Wilson loop on 
the lattice. Throughout, we only consider /). The Creutz ratio is interpolated between 
values of I to define a related function x{R) with R = R/a continuous. The interpolation 
requires a matching of values at the points where R is half-integral: 

.An .a 1 In , W{R+^,R+l)W{R-l,R-l) ^ ^ 

X(R) = x{R+ -,R+ -) = - ^^ — - — r-^ 2.2 

the logic being that the Wilson loops appearing in the Creutz ratio have average size R. In 
fact, taking the classical continuum limit, one finds that the Creutz ratio is a finite difference 
approximation : 



R'xiR) ^ -RT^^\nW,,MT) 



(2.3) 

T=R 



where Wcont.(-R, T) ^ W{R, T) is the Wilson loop in the continuum language. T here should 
not be confused the temporal extent of the lattice. 

In the Creutz ratio method for determining g{L), we choose 

r = R/L = fixed. (2.4) 

Note also that we have to deal with a fermion mass nig, which we take to be the PCAC 
mass. Our approach will be to measure g{L) at nonzero and then extrapolate to the 
ruq = limit. 

The important thing that Bilgici et al. have noted is that at one loop in lattice pertur- 
bation theory, including the effect of bosonic zero modes, 

^'xi-ioop(^, a\R) = k{r, N)gl N = L/a (2.5) 



1 In our analysis we take temporal extent T = 2L, to allow for accurate measurement of the partially 
conserved axial current (PCAC) mass and reuse of the lattice configurations for studies of the spectrum. 
In fact, this is an advantage of the Creutz ratio method over Schrodinger functional studies. 
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which defines the quantity k{r,N). Here, go is the bare couphng. It was found in 
that in the hmit of large N, k{r, N) — )■ k{r), so that it is an L and a independent quantity 
in the continuum hmit. Thus we write k{r) in what follows, and will define the running 
coupling in terms of this — ?■ oo normalization factor. Fermions enter x at two-loop order, 
so k does not depend on niq, since it is defined by the one-loop expression. 

Thus we can make an L-dependent nonperturbative definition of the running coupling g 
using the value of x{R) obtained from a simulation with bare coupling go:"^ 

fiL) = lim -^R^xiL, a\R), Ra/L = R/N = r = fixed (2.6) 

a^o k[r) 

Note that as we shrink a, in order to hold L fixed the number of lattice sites N increases, 
and so must the size of the Wilson loops R. Since these quantities fall off exponentially 
with increasing R (relative to noise) and larger lattices are significantly more expensive to 
simulate, this becomes a demanding computation, particularly with dynamical fermions. 
To get a handle on g'^{L), we can study the scaling of the function: 

(2.7) 



g^(N,r,L)^-^R'xiN,a\R) 



R=rN 

On the right-hand side, the nonperturbative xi^^^^l^) is evaluated with system size L. 
Clearly if L is held fixed while N is increased, the lattice spacing is decreasing toward 
the continuum limit. Thus g'^{N,r,L) of (12. 7p defines the running coupling in a particular 
(nonperturbative) scheme: 

g\r\L)= hm g\N,r,L). (2.8) 

Af— s>oo 

The constant r is part of the renormalization scheme. 

A step-scaling analysis is used to follow the running of the coupling g with the scale L. 
We will not actually use step-scaling in order to search for the presence of an infrared fixed 
point. Rather, there is a simpler, qualitative behavior that we are looking for, as we now 
describe. 

B. Expectations 

In a confining theory, g'^{xN,r,xL) > g'^{N,r,L) for x > 1. This type of relationship 
means that if we hold /3 fixed, the measured value of g'^{N,r,L) increases as N increases, 



^ Often below we will use the lattice coupling /3 = 4/(7o- 
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since L = Na will increase proportionately. This leads to a behavior sketched in Fig. [T] In 
fact, it is easy to convince oneself that a step-scaling analysis applied using the curves in 
Fig. [T]will lead to a g^{L) that increases with L. 

On the other hand, if there is an infrared fixed point, the ordering of curves with in- 
creasing reverses once the bare coupling (3 goes past the fixed point (/3 < /3*(A^)), and 
g'^{xN,r,xL) < g'^{N,r,L) occurs for x > 1. Here we note that the fixed point /3* depends 
on A^, due to lattice artifacts. Indeed, such a "backwards running" was seen at strong 
couplings in the SU(3) theory with sixteen fundamental flavors in the work 20|, |2l(|. We 
expect something like this to also occur in the present theory at strong coupling, under the 
assumption that there is an infrared fixed point. Thus we expect to obtain a behavior like 
the one sketched in Fig. [21 Note that the curves do not all cross at the same point, which 
is why /3* is a function of A^. Notice that there is a region of small /3 where the ordering of 
curves has clearly reversed, and that if we did a step scaling analysis in this region, the flow 
would be backwards. 

From the vantage point of the two sketches that we have just presented, it is clear that a 
full blown step scaling analysis is not needed in order to detect the presence of an infrared 
fixed point. All that is needed is a plot of g'^{N, r, L) versus /3 for the various A^. In order to 
decide whether or not there is an infrared fixed point, the question that we need to answer 
is just this: "Which of the two figures does the plotted data look like. Fig. [1] or Fig. [2]?" 
In fact even less is needed: it is already known from perturbation theory at weak coupling 
that the curves are ordered with respect to L as shown at large (3. So what we need is just 
to show the reversal of ordering with respect to L at small /3. This will be the focus of our 
analysis. 



III. ANALYSIS 



A. Smearing 



We follow Bilgici et al. and apply APE smearing 2J, |25| for the links that are used in 
the calculation of Wilson loops: 



Uf,{x) P{U^{x) + a '^[Uu{x)U^{x + aO)Ul{x + afi) 
+Ul{x — aO)Un{x — ai')Uy{x — aO + afi)]} 



(3.1) 




where P projects to a unitary matrix, which for the SU(2) group that we work with is 
accomphshed simply by 

A ^ ^ (3.2) 

for a matrix A. We take the smearing parameter to be a = 0.5. We find that two smearing 
steps works well for yielding the perturbatively required increase in g'^{N,r,L) with N at 
weak bare coupling. For this reason we use two smearing steps throughout this work. 



B. Interpolation 

In order to obtain a value of for r = 1/4 with arbitrary N, it is necessary to interpolate 
to values of R that are not half-integral. For this purpose Bilgici et al. introduce a quadratic 

hypothesis: 

kg'^{R) = R?x{Ro) Co + ciR + C2R? (3.3) 

However, the obtained value is somewhat different if a cubic hypothesis is chosen. Further- 
more, the result depends on which points are included in the fit. Here the considerations 
are that the smallest values of R are degraded by the smearing, and the largest values of R 
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start to reflect periodicity. We have adopted the strategy of varying all these choices and 
then using the variation of kg'^ as a measure of systematic error in the method. The mean 
value of kg"^ is used as our final estimate. Further details are given in Appendix O 



C. Results 

Now that we have outlined how we obtained the values of k{r)g'^{N,r, L), we proceed 
to discuss our results.^ We have measured kg"^ for lattices N = 10, 12, 16, 20, 24 with bare 
coupling P = 2.25. Note that this value of /3 is at weaker coupling than the one where a bulk 



phase transition occurs, /3 ~ 2, as was found in 26|]. It is also at a stronger coupling than 



the P where there was some evidence for fixed point behavior in the Schrodinger functional 



study 151 . Thus this point is continuously connected to the continuum limit through finite 
mass, but is expected to be on the strong side of the infrared fixed point. That is, we 
have reason to suspect backwards running for the choice /3 = 2.25. We have performed our 



^ We do not need the actual value of k{r = 1/4) since it is a fixed constant independent of N. However, in 
Appendix |B] we consider an alternate definition of the running coupling where k depends on TV, and in 
that case we do determine and divide by it. 
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simulations for five values of the bare mass moa, shown in Table [H We have measured the 
PCAC mass m^a from the largest lattices (A^ = 24), where there is the least systematic error 
from finite volume. These are also given in Table HI We then extrapolate k{r)g'^{N,r, L) to 
the rriq = limit with quadratic and cubit fits to the data. In the case of the quadratic fit, 
we also consider the case where only the lightest four masses are included in the fit. The 
results for the measured values of kg"^ are summarized in Table [TTl and the results of the zero 
mass fits are presented in Table lllli The results of Table [TTl are displayed along with the 
quadratic fits to all data in Fig. [31 

At large mass one sees the behavior characteristic of asymptotic freedom: 

g\24) > /(20) > ^7^(16) > g\l2) > /(lO). (3.4) 

In the case of the quadratic fit to all data, when we extrapolate to the chiral limit the trend 
reverses: 

/(24) ^ /(20) < /(16) < /(12) ^ /(lO). (3.5) 

Thus for this extrapolation we find backward running in the massless limit. However, in 
the other two extrapolations of Table IIII| what we see is that there is no clear pattern, but 
instead a rough equality: 

^7^(24) ~ /(20) ~ /(16) ~ /(12) ~ /(lO). (3.6) 

Thus what one has in this case is evidence for an absence of running. 

We note that the P^r degree of freedom (d.o.f.) is quite large for most of the fits. The 
error estimates take this into account but the implication is that the mass dependence of kg"^ 
is complicated and not a simple polynomial. In such a case on can question the reliability of 
the zero mass extrapolations. Appeal must finally be made to the results at rriQa = —1.18, 
which shows a rough equality between the values of kg"^. We expect that this will also hold 
for the massless limit. In conclusion, our results seem to favor an absence of running, though 
there is a hint of backwards running from one of the extrapolations. 

The behavior that we have observed is consistent with the existence of an infrared fixed 
point. We regard our results as suggestive that Minimal Walking Technicolor does not 
actually walk, but is instead inside the conformal window. This is supportive of the findings 
of Schrodinger functional studies, but now by a different method. 
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rriQa 


rriqa 


-1.000 


0.448406(5) 


-1.100 


0.236337(4) 


-1.165 


0.090917(2) 


-1.175 


0.066873(6) 


-1.180 


0.054890(3) 



TABLE I: PCAC mass obtained from 24^ x 48 lattices. 




0.1 0.2 0.3 0.4 0.5 



rriqa 

FIG. 3: Values of kg^ for /3 = 2.25, r = 1/4. At large mass, increases with increasing so 
that the curves are ordered top to bottom, N = 24, 20, 16, 12, 10. However, in the massless limit, 
the ordering of the curves shows a trend of reversal, indicating backward running. 

IV. CONCLUSIONS 

We have found that the qualitative behavior of the couphng g'^{N, r, L), defined through 
the Creutz ratio, is consistent with the existence of an infrared fixed point. We arrived at 
this conclusion by arguing that two different behaviors are possible, sketched in Figs. [Hand 
[2J We compared our data for g'^{N, r, L) to these figures and find that as the fermion mass 
is reduced toward zero, the ordering of the points either reverses at strong bare couplings or 
is roughly constant, corresponding to the region where there is a crossing of curves in Fig. |2l 
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moa = -1.000 


moa = —1.100 


moa = —1.165 


moa = —1.175 


moa = —1.180 


10^ X 20 


0.6448(67) (0) 


0.4832(18)(0) 


0.38394(48) (0) 


0.36808(87) (0) 


0.3630(10) (0) 


12^ X 24 


0.7604(6) (32) 


0.5154(2)(32) 


0.3884(2)(16) 


0.37183(16)(66) 


0.35661(19)(35) 


16^ X 32 


0.9820(10)(37) 


0.6114(8)(31) 


0.3954(3)(18) 


0.3647(2)(13) 


0.35906(20) (97) 


20^ X 40 


1.2621(12)(21) 


0.6998(11)(28) 


0.4018(3)(13) 


0.3703(4)(22) 


0.3474(3)(17) 


24^ X 48 


1.6510(15)(7) 


0.8109(15)(12) 


0.4192(5)(10) 


0.3870(24) (7) 


0.3623(8)(18) 



TABLE II: Results for kg^ for the various sizes of lattices and bare masses. The first error in kg^ 
is statistical, determined by jackknife analysis of fits to the Crcutz ratios as described in the text. 
The second error in kg^ is systematic, obtained from varying the method of interpolation used to 
obtain the value at R = 0.25(L/a). This is described in more detail in Appendix blah. 





Quadratic, last 4 pts. 


Quadratic, all 5 pts. 


Cubic, all 5 pts. 




kg' 


x'/d.o.f. 


kg' 


x'/d.o.f. 


kg' 


x'/d.o.f. 


lO'^ X 20 


0.3312(49) 


2.52 


0.3288(20) 


1.75 


0.3321(58) 


2.44 


12^ X 24 


0.297(17) 


23.6 


0.3130(95) 


25.7 


0.291(20) 


21.3 


16^ X 32 


0.319(20) 


12.7 


0.295(11) 


19.0 


0.327(23) 


11.4 


20^ X 40 


0.284(15) 


5.04 


0.2651(85) 


7.88 


0.288(18) 


5.36 


24'^ X 48 


0.313(18) 


7.12 


0.268(11) 


29.5 


0.321(22) 


7.81 



TABLE III: Results for the zero mass extrapolation of kg' for the various sizes of lattices. 

While it would be disappointing from a phenomenological point of view if Minimal Walk- 
ing Technicolor is actually conformal in the infrared, it is nevertheless significant that by 
more than one method lattice field theory technique evidence has been obtained for the ex- 
istence of an infrared fixed point. Furthermore, conformal theories are inherently interesting 
and more detailed studies of the present gauge theory are worth pursuing on the lattice. 

Two refinements to the present study could be performed, which would hopefully yield 
firmer conclusions. First, smaller PCAC masses should be simulated so that a stronger case 
for what happens in the chiral limit can be made. Second, larger volumes such as N — 32 
should be simulated, so that the trends with respect to N can be amplified, hopefully 



10 



clarifying whether or not backwards running actually occurs. We are presently pursuing 
these two improvements and will report on them in the future. 
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Appendix A: Computing methods and resources 

Simulations are performed using Hybrid Monte Carlo [27[, which is the standard method 
for dynamical fermions. Our lattice action is Wilson fermions and a plaquette gauge action. 
Our sampling of configurations is separated by five molecular dynamics time units, and we 
use 1000 samples on all lattices excepting = 24, niQa = —1.175,-1.180 where we used 
700 samples in each case. (This was the most expensive part of the computation.) The 
parallel simulation and analysis code that we use is a modification of the Columbia Physics 
System, in order to incorporate SU(2) gauge group and adjoint fermions. It was originally 



developed for the study of A/" = 1 super- Yang-Mil 



but has recently been used for 



the Minimal Walking Technicolor gauge theory [26|, |29[. All simulations were performed on 
IBM BlueGene/L computers at RPI, both an on-campus rack and one to four of the sixteen 
racks at the Computational Center for Nanotechnology Innovation. Our larger lattices were 
simulated using an entire rack, while our smaller lattices used smaller partitions. Suffice it 
to say, significant computing resources were dedicated to this project over a period of about 
one year. We estimate that 7 million BlueGene/L core hours were expended for the research 
summarized in this article. Smearing to obtain the Wilson loops that are used for our Creutz 
ratios was also performed in parallel, though this analysis step consisted of only a small 
fraction of our computing time. Finally, we computed the correlation functions necessary 
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for obtaining the PCAC mass using the parallel resources, which again only consumed a 
small fraction of the computing time. 

The remainder of the numerical analysis for this project was performed on a single work- 
station. This involved averaging the data with a jackknife analysis of fits in order to obtain 
errors in our values of kg"^. Similarly, the final steps of the analysis to obtain PCAC masses 
were performed with jackknife error estimation. 



Appendix B: Alternative definition 



Since there is no unique definition of the running coupling at finite A^, one could also use 
the value of k at finite A^. That is, 

1 



g'{N,r,L) 



R'x{L,a\R) 



(Bl) 



k{r, N)' 

For this, since we choose r = 1/4 it is necessary to interpolate k{r,N) between values of r 
that are allowed for a given A^ (we use a quadratic fit at fixed A^). For the determination 
of k, we use the results of 30|, which leads to the following formulae. First their is the 
contribution of Wilson loops which do not have links in the temporal direction: 

2 



cf\l,J) 



m {4v^\p(pW 



- II 



J-1 
A;=0 



+ (/ ^ J, 1 O 2; 



(B2) 



Nc is the number of colors; A"c = 2 in our case. The quantities |P(p)p and V are defined 
below. A similar formula holds for the Wilson loops with two links in the temporal direction: 

2 



cf)(J,J) 



N',-l{ 1 



AV ^ I P{-n 



Ar2 l4vA^|p(p)|2 



JPol 
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J-1 



J2 e"^'"' 



k=0 



(B3) 



In these expressions, the momenta take the following values, given that there are A^ sites in 
the spatial directions and 2A^ sites in the temporal direction {i = 1, 2, 3): 



Po = non/N, p, = 2ni7i/N, no = 0, . . . ,2N - 1, ni = 0,...,N -1 



(B4) 



Also, we have the definitions: 



(B5) 



12 



0.15 



0.10 — 



0.05 — 



0.00 



X N=10 

o N=12 

□ N=16 

* N = 20 

« N = 24 



0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



FIG. 4: k{r, N) for the lattices that we study. 



The zero-model contribution must also be included: 



cf) =ci(p = 0) + ci' 



(0) J2) 



ci{p = 0) + cf\ ci(p = 0;/,J) 



12m V 



Finally, the Creutz ratio at this order of lattice perturbation theory is: 



1, 



l-loop 



glN,;-[cf\l, I) + cni, I) + cY>{I - 1, / - 1 



.(0) 



(2), 



+cf\l - 1,/ - 1) - c\'>{I,I- 1) - cY>{I -1,1)- 2cr(/,/ - 1)] (B7) 



.(2) 



(2), 



.(0), 



(B6) 



where Qq is the bare coupling. From this we extract k at finite N 

2 



A;(r, iV) = ^ (^/ - i ) x(/,/)i-ioop, r 



I-l]/N. 



(B8) 



In Fig. m we show k{r, N) for the lattices that we study, taking the range 1 = 2,..., N/2. We 
note that our values of k{r,N) are in agreement with those of 18|, ll9|, once the difference 
in color factors (A^^ — 1)/Nc is taken into account (they have A^^^ = 3). 

We have calculated by this alternate approach, for the case of quadratic extrapolation 
fitted to all data points, and the results are summarized in Table IIVI The behavior of 
decreasing coupling with increasing N seems to persist with this modified definition of the 
running coupling. 
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0.111952 


2.937(23) 


0.108871 


2.875(87) 


0.103207 


2.85(10) 


0.0995411 


2.663(85) 


0.0971746 


2.76(14) 



TABLE IV: calculated by the method where different values of k are used for each N . 



Lattice 


-^start 




Fit 


10^ X 20 


2 


5 


2 


12^ X 24 


3 


7,8 


2,3 


16^ X 32 


3,4 


7,8 


2,3 


20^ X 40 


3,4,5 


7,8 


2,3 


24^ X 48 


4,5 


9,10 


2,3 



TABLE V: Choices of parameters that were varied for the fit that is used in the interpolation to 
obtain kg'^{R = N/4.). 

Appendix C: Systematic error estimates 

As described in the main text, an interpolation of kg'^{R) to the value R = N/A is 
required, since x{^^ ^) only takes values at integer I and R = I — j. This allows for various 
choices for the range of the fit, /start, • • • , hnd- One can also choose to fit with a quadratic 
(Fit=2) or cubic (Fit=3) hypothesis for kg'^{R), as a function of R. We have varied these 
options, constrained by what is possible on the various lattices, and have used the standard 
deviation of the interpolated values as a measure of systematic error. Because of smearing 
it is important that we choose /start as large as possible on a given lattice. In Table |V] we 
present the various parameter choices that were used. The mean of the interpolated values 
is used as our estimate for kg"^. 
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